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SECOND-ORDER FLUCTUATIONS AND CURRENT ACROSS 
CHARACTERISTIC FOR A ONE-DIMENSIONAL GROWTH 
MODEL OF INDEPENDENT RANDOM WALKS^ 

By Timo Seppalainen 

University of Wisconsin 

Fluctuations from a hydrodynamic limit of a one-dimensional 
asymmetric system come at two levels. On the central limit scale 
n}^^ one sees initial fluctuations transported along characteristics 
and no dynamical noise. The second order of fluctuations comes from 
the particle current across the characteristic. For a system made up 
of independent random walks we show that the second-order fluc- 
tuations appear at scale n^^* and converge to a certain self-similar 
Gaussian process. If the system is in equilibrium, this limiting process 
specializes to fractional Brownian motion with Hurst parameter 1/4. 
This contrasts with asymmetric exclusion and Hammersley's process 
whose second-order fluctuations appear at scale n^^^, as has been 
discovered through related combinatorial growth models. 

1. Introduction. An interface model defined in terms of a height function 
on an integer lattice is a stochastic process at = {o't{x) :x € Z*^}, where the 
value at{x) is interpreted as the height of the interface over site x. The 
random variables ctj (x) move up and down according to random rates whose 
momentary values depend on the height values in some neighborhood around 
site X. A hydrodynamic scaling limit is a type of law of large numbers for 
these systems. The conventional statement is that under a suitable scaling of 
space and time, the entire space-time random evolution {crt{x) : x € Z'^, t>0} 
converges to the solution of a differential equation. 

When the system is asymmetric, in the sense that there is an average 
drift either up or down, the typical law of large numbers for the system is 
of the following form. The result is for a sequence of processes a", where 
n = 1, 2, 3, . . . is the scaling parameter. The statement is that as n ^ oo, the 
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random position n ^a^^{[nx]) converges to a function u{x,t) that solves a 
Hamilton-Jacobi equation 

(1.1) ut + f(yu) = o. 

In the macroscopic description, / gives the local velocity of the height as a 
function of the local gradient. A necessary hypothesis for this type of law 
of large numbers is that as n — > oo , the sequence of scaled initial states 
n~^(To([ny]) in some sense converges to a function UQ{y). The function uq 
then serves as the initial data for the equation ut + f{Vu) = 0. The equation 
and uq uniquely determine u{x,t) at all later times t. Depending on the 
situation, there might be additional assumptions on the distributions of a^. 
Examples of results for various models can be found in [12, 14, 16, 17, 19]. 
For general accounts of hydrodynamic limits, we refer to the monographs 
[11] and [20], and to the lectures [4] and [21]. 

The fluctuation question for asymmetric systems has so far found an- 
swers only in the one-dimensional situation, where the following picture has 
emerged. Suppose the initial conditions satisfy a central limit theorem of 
the type 

*Co(y), yGR, 



n 

with a continuous limiting process Co- This situation arises naturally when 
the initial increments (Tq{x) — Uq (x — 1) are independent with slowly varying 
bounded means and variances. Then at later times a weak limit 
(1.2) <,([nx])-nn(:r,t)^^ 



holds. The process Ci^it) is a deterministic function of the initial process 
Co- More specifically, the value Ci^^t) is determined by the values Coiy)^ 
such that a generalized characteristic of the p.d.e. (1.1) emanating at (y,0) 
reaches {x,t). Qualitatively, a crucial feature is that there is no dynamical 
noise visible at this scale n^^^ . These types of results have been proved for the 
exclusion process under various hypotheses [7, 8, 13] and for Hammersley's 
process [18]. 

The motivation of the present paper is to describe fluctuations that lie 
"beyond" the trivial fluctuations transported by the characteristics that ap- 
pear in (1.2). This second level of fluctuations appears when the first-order 
fluctuations of (1.2) are suitably subtracted off or when the initial conditions 
are deterministic. In the asymmetric exclusion and Hammersley settings, one 
should find some kind of fluctuations at the n^/^ scale. The results of Baik- 
Deift-Johansson [2] and Johansson [10] can be interpreted as fluctuation 
results for Hammersley's process and the exclusion process with special de- 
terministic initial configurations. So for exclusion and Hammersley's process. 
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the aim would be to generalize those v}^'^ results to other initial conditions. 
What is special about exclusion and Hammersley's process is that these 
"second-order" fluctuations arise through natural growth models that are 
amenable to the powerful combinatorial and analytic machinery of [2] and 
its descendants. 

In the present paper we start another direction, the investigation of these 
phenomena in other models besides exclusion and Hammersley's process. 
Two questions arise. Do the trivial n^/^ fluctuations transported by char- 
acteristics appear in other asymmetric models? What then would be the 
second-order fluctuations, especially if there is no combinatorial growth 
model present that would lead to the n^/^ fluctuations and the random 
matrix connections? 

A simple model is one where the increments of the height function come 
from independent random walks. This case we analyze in the present paper. 
The initial increments are taken independently with slowly varying means 
and variances. Exactly as for exclusion and Hammer sley, on the central limit 
scale n^/^ we have the initial fluctuations transported by characteristics. 
Then we find the next order of fluctuations on the scale v}/^. In the limit 
these fluctuations are described by a certain Gaussian process Z{y,t), where 
(y,0) is the initial point of the (now unique) characteristic that reaches 
{x,t). The covariance of Z{y,-) is determined by the mean and variance of 
the initial increments around the macroscopic point y. Imprecisely speaking, 
the random height expands as 

<i (M) ~ nu{x, t) + n^/'^Q{x, t) + n^'^Z{y, t). 

The processes Z(y, •) are independent for distinct initial points y. 

In the special case when the height increments are in equilibrium, for a 
fixed y, the process Z{t) = Z{y,t) specializes to fractional Brownian motion 
with Hurst parameter j. These second-order fluctuations turn out to be 
the same as the particle current across a characteristic of the macroscopic 
equation. Hence, the juxtaposition in the title of the paper. 

In the next section we describe the random walk model, its hydrodynamic 
limit and the two levels of fluctuations. For the sake of comparison, we 
include a brief section on the fluctuation picture of Hammersley's process. 
We show that for Hammersley's process the second-order fluctuations are, 
at most, of order n^^^ logn. This bound is valid also for shock locations. The 
last two sections prove the theorems. In the proofs, C, Co, Ci,... denote 
constants whose actual values may change from line to line. The set of 
natural numbers is N = {1, 2,3,...}. 

2. The random walk model. We consider a model of an interface whose 
height differences between neighboring sites are defined by independent ran- 
dom walks on Z. The state of the system at time t is a height function 
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(7f : Z — > Z. The height function is nondecreasing in space, so the increments 

(2.1) T]t{x) = at{x) - at{x - 1) 

are nonnegative integers. 

The randomly evolving height function is constructed as follows. Let 
{Xi{t) : i G X} be a countable collection of independent continuous-time ran- 
dom walks on Z. The jump rates of the random walks are given by a prob- 
ability kernel {p{x) :x G Z}. In other words, the assumption on p{x) is 

Y.p{x) = i 

and the common transition probability of the random walks is 

oo —tj.k 

P[X,{s + t)= y\X,is) =x]= ptix, y) = Y^ ^P^'Kv " 

k=o ^■ 

where 

pW(z)= p{xi)p{x2)---p{xk) 

Xl+X2-\ \-Xk=Z 

is the /c-fold convolution of the kernel p{x). 

Given an initial height function ctq = {(To(x) : x G Z} defined on some prob- 
ability space, define the initial increments r]o{x) = (7o{x) — ao{x — l). Choose 
the initial positions of the random walks so that site x contains r]Q{x) par- 
ticles: 

(2.2) J2^{X,iO) = x} = mix). 

Once the initial points {Xj(0) :i G X} have been specified, define the subse- 
quent evolutions {Xi{t) — Xj(0) : i G X, t > 0} as an i.i.d. collection of random 
walks on this same probability space, independent of a^. Define the current 
Jt{x) as the (net) number of particles that have moved from (— oo,x] to 
[x + l,oo) in time interval [0,t], 

Jt{x) = ^ i{x,(o) < X < Xi{t)} - HMt) <^< MO)}- 

i i 

The height function at time t is then defined by 

(2.3) at{x) = ao{x) - Jt{x). 

In other words, the interface height at x moves down one step with every 
particle that jumps from (— oo, x] to [x + l,oo), and up one step with every 
particle that does the opposite. The increment variables rit{x) defined by 



CURRENT ACROSS CHARACTERISTIC 



5 



(2.1) also serve as the occupation variables of the random walks: from (2.1)- 

(2.3) one can derive 

(2.4) r^t{x)=J2l{X,{t)=x}. 

i 

We can describe the evolution of at directly in terms of the rates, without 
reference to the random walks. Given a height function a (a nondecreasing 
function Z — > Z), and x,£ £ Z, define a new height function cr^'^ by 

a{y)-l, x<y<x + £-l, if^>o- 
_ J lcr(y), if y<x or y>x + i, - ' 



(2.5) a^'^y) 



a{y) + l, x + i<y<x-l, if£<o 
cr(y), if y < X + £ OT y > X, 



The dynamical rule for the height process is this: if the current state is a, 
then for each a; G Z and i G Z, at rate p{£){a{x) — a{x — 1)), the process 
jumps from a to <t^'^. If £ = 0, there is actually no change: o"^'" = a. 

Now some assumptions. We give them in three groups, first the assump- 
tion on the kernel p{x) and then the assumptions on the sequence of initial 
height functions {ctq } that determine the hydrodynamic limit setting. 

Assumption A. For the random walk kernel, we assume that, for some 
6>0, 

(2.6) e^'^Pix) < oo for 1^1 < 6. 

The purpose of Assumption A is to enable us to use standard large devi- 
ation bounds on the random walks. 

Assume given a sequence of initial height functions ctq , random or de- 
terministic, defined on some probability space. The positive integer param- 
eter n will tend to oo in the results. Define initial occupation variables 
rjl^ix) = a^{x) - aJJ^{x - I). 

Assumption B. Assume that for some nondecreasing function uq 
on R, and all y G R, 

(2.7) lim an {[ny]) = uo(y) in probability. 

n — ^oo 

Assumption B is for the hydrodynamic limit. For the fluctuation results 
we need stricter control of the initial conditions, as in the next assumption. 

Assumption C. For each n, the initial occupation variables {f?o (x) : x £ 
Z} are independent, with a uniformly bounded sixth moment: 

(2.8) sup E[r]^{xf] < oo. 

nSNjXeZ 
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Let 

p]^{x) = Er]^{x) and = Var[7?^(x)] 

be the mean and variance of the initial occupation variable rj^^x), x S Z. 
Let uq be the function specified in Assumption B, set pQ = u'q^ and let vq 
be another given nonnegative function on R. Assume both po ^ind vq are 
bounded. The means Pq{x) and variances Vq{x) approximate the functions 
po and vo in the following precise sense: 

For each y E R, there exist positive integers L = L{n) such that n~^^^L{n) - 
and for any finite constant A, 



(2.9) 



lim 

n— >oo , 



sup 



n 



1/4 



L(n) 



L(n,) 

E 



Poib^'y] + + j) - po{y) 



0. 



\m\<A-^nlog n 

The same assumption holds when Pq and po ^ire replaced by and vq. 

Throughout the paper, L = L{n) denotes the quantity specified in the 
assumption above. The awkwardly complicated assumption (2.9) is made 
to accommodate both random and deterministic initial conditions. Given a 
function po, the expectation of a random 77o(x) can, of course, agree exactly 
with Po(f )i but a deterministic 7?o (^) cannot unless Po{^) is integer- valued. 
Here are two basic examples of initial conditions that satisfy Assumptions 
B and C. 



Example 2.1 (Random initial conditions). The functions po = Uq and 
"^o are bounded, nonnegative and satisfy a local Holder property: for each 
bounded interval [a, 6], there exist (3 = P{a,b) > 1/2 and C = C{a,b) < oo 
such that 

(2.10) \po{x)- po{y)\ + \vo{x)-vo{y)\<C\x-y\f^ for all a;,y G [a,6]. 

For each n, let {'^^{x) :x G Z} be independent, satisfy assumption (2.8) and 
have 



(2.11) 



n 



En^{x)=po{-] and VarK(x)] 



X 

vo\- 
n 



Additionally, the variables {ctq (0)} are chosen so that (2.7) holds for y = 0. 
Then (2.7) is satisfied for all y G R. 

Example 2.2 {Deterministic initial conditions), po = Mq is bounded, 
nonnegative and satisfies the Holder condition (2.10) with (3 = (3{a, b) > 1/2, 
and Vq is identically zero. Define deterministic initial occupation variables 
by 

m — 1' 



(2.12) 



nuQ 



m 



n 



nuo 



n 
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where, as throughout the paper, [x] = inax{k G Z : /c < x} denotes the inte- 
ger part of a real x. The function uq is nondecreasing so each r]Q{m) is a 
nonnegative integer. And finaUy, the variables {cq (0)} can be random or 
deterministic, but must satisfy (2.7) for y = 0. 

We set the stage with the hydrodynamic limit and the first-order fluctu- 
ations. The first two moments of the random walk kernel appear in various 
parts of the results. We denote these by 

b = ^'^xp{x) and K2 = x^p(x). 

X X 

Let 

(2.13) u{x,t) = uo{x -bt). 
It is the solution of the linear transport equation 

(2.14) Uf-F^ii^,. = 0, u(x,0) = uo(x). 

Theorem 2.1. Assume Assumptions A-C. Then, for each {x,t) £ R x 
[0,oo), 



(2.15) limn a^^^dnx]) = u{x,t) in probability. 
Given {x,t), lety = x — bt. Then 

(2.16) lim l<t{[nx])-nu{x,t) _ a^{[ny]) -nupm ^ ^ pMUty. 



Furthermore, assume the specific situation of Example 2.1 and the nor- 
malization uo(0) = o'q{0) = for all n. Then it is possible to construct the 
processes on the same probability space with a two-sided standard Brow- 
nian motion B{-) such that these limits hold in probability, for all {x,t), 
y = X — bt: 

(2.17) lim <.(M)-n^(^.o^^(-r,^(,),,y 



A two-sided Brownian motion B{-) is constructed by taking two indepen- 
dent standard Brownian motions Bi and B2 on [0, 00) and setting 

Bi{t), t>0, 
-B2i-t), t<0. 



B{t) 



The integral Jq vo{s) ds in (2.17) has to be interpreted with a sign, in other 
words, 

y j-o 
VQ{s)ds = — Vo{s)ds for y < 0. 
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The characteristics of (2.14) are straight hnes with slope b, so x = y + bt is 
the characteristic starting at (y,0). Limits (2.16) and (2.17) say that on the 
central limit scale n^/^ , fluctuations from the hydrodynamic limit consist of 
initial fluctuations rigidly transported along the characteristics, without any 
contribution from dynamical noise. 

The purpose of this paper is to describe the "second-order" fluctuations 
that appear beyond the trivial fluctuations of Theorem 2.1. Fix y G R. Let 

(2.18) Yn{t) = a^ny] + [nbt]) - a^{[ny]). 
Since u{y + bt,t) = UQ{y), also 

(2.19) <i([ny] + [nbt]) - nu{y + bt, t) = a^{[ny]) - nuo{y) + y„(t). 

So Yn{t) represents the difference between the fluctuation experienced by 
the process at space-time point {y + bt, t) and the fluctuation at the initial 
point (y,0) of the characteristic. n~^/^l^,(t) is exactly the difference that 
appeared in (2.16). 

There is another way to look at Yn{t), directly in terms of the particles. 
Write J" for the current of process a" and X^(t) for the random walks in 
the construction of erf: 

Y^{t)=a^,{[ny] + [nbt])-a^{[ny]) 

= a^iiny] + [nbt]) - j;i,{[ny] + [nbt]) - a^i[ny]) 

[ny] + [nbt] 

= E V^{m)-JMny] + [nbt]). 

m=[ny]+l 

Switching back to the random walks and cancelling gives 

Ynit) = E HKiO) > [ny] + l,X^int) < [ny] + [nbt]} 

(2.20) 

- E l{^r(0) < [ny],Xnnt) > [ny] + [nbt]}. 

i 

This counts the net number of particles that have moved from the right 
side of the characteristic to the left side during time interval [0,nt]. In other 
words, Yn{t) also represents the negative of the current across the charac- 
teristic. 

Our main theorem is the distributional limit of Y^. Assumption B is not 
relevant for the limit of Y^. The previous paragraph already showed that 
even though Yn{t) was defined in terms of the height functions in (2.18), it 
is actually determined by the increment process 77" . We included Assumption 
B in the earlier discussion only to give the complete hydrodynamic picture. 
Also, the approximation assumption (2.9) is needed only for the particular 
y that appears in the definition of Yn- 
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Theorem 2.2. Fix y G R and define Yn{t) as in (2.18). Assume As- 
sumptions A and C. Then as n ^ oo, the process n"^/^l^(-) converges 
weakly on the space L'r[0,oo) to the mean-zero Gaussian process Z{-) with 
covariance 



EZ{s)Z{t) = . -±{po{y){./7Tt - VsVt-sAt) 

V ZTT 

(2.21) ^ 

+ vo{y)i^/^+Vt- Vl+t)}. 

For the increment process rjt = {r]t{m) -.rnGZ}, i.i.d. Poisson distributions 
are equilibrium distributions. If the mean of the Poisson is p, then po{x) = 
vq{x) = p for all X gH. The covariance in (2.21) then simplifies to 

(2.22) EZ{s)Z{t) = pJ!^{^/^ + Vi- Vt^) fors<t. 

V 27r 



This is the covariance of fractional Brownian motion with Hurst parameter 
H = J, normalized by EZ{1)'^ = p^J^K/^pK. Standard fractional Brownian 
motion would have EZ{1)^ = 1. Let us state this special case as a corollary. 



Corollary 2.1. Letr]t be an equilibrium process whose occupation vari- 
ables {r]t{x) :x G Z} are i.i.d. Poisson with mean p. Let ry" = rjt for each n. 
Then ra-i/4y„(-) converges weakly on the space D^[0,oo) to fractional Brow- 
nian motion with covariance given in (2.22). 

Another special case worth highlighting is that of deterministic initial 
height functions described in Example 2.2. In that case the limit in (2.17) 
is zero because the second fraction in (2.16) vanishes in the limit. (2.18) 
gives the actual fluctuations from the hydrodynamic limit because aQ{[ny]) 
is deterministic. We also omit the short derivation of this corollary from 
Theorem 2.2. 

Corollary 2.2. Let uq be a nondecreasing function on R with 
bounded derivative po. Assume pQ satisfies the local Holder condition (2.10). 
Define deterministic initial height functions by ctq (?7i) = [nuo( — )], m G Z. 
The function vq is now zero. Fix y E R. Then the process 

<ti[ny] + [nbt])-nu{y + bt,t) 



n 



1/4 



:t >0 



converges weakly on Dr.[0,oo) to the mean zero Gaussian process Z{-) with 
covariance 

EZ{s)Z{t)=pQ{y).[^{^/7Tt-Vt^s) fors<t. 

V ZTT 
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Fractional Brownian motion has stationary increments, but the general 
process Z{t) with covariance (2.21) does not unless poiv) =vo{y). One can 
check that for a fixed h> 0, E[{Z{t + h) — Z{t))'^] is strictly decreasing with 
t if uo(y) > Poiy), and strictly increasing if uo(y) < Po{y)- A bound 

E[{Z{t)-Z{s)f]<C{t-s)'/' 

is valid for all < s < t. Since the increment is mean-zero Gaussian, it follows 
from Kolmogorov's criterion that the process Z has continuous paths. 
The process Z is self-similar with index |, which means that {Z{at) -.t > 

0} = {a}^^Z{t) :t > 0}, as is immediate from the form of the covariance. 

Next we address the joint distribution of processes Yn{-) from several 
initial points y. Write Yn{y,t) for the random variable defined by (2.18) or, 
equivalently, by (2.20) to display its dependence on y. Write Z{y,t) for the 
Gaussian process with covariance given in (2.21). 

Theorem 2.3. Assume Assumptions A and C. Let yi < y2 < ' ' ' < Vk be 
points on R. Then as n ^ oo, the joint process n~^/^{Yn{yi, ■),Yn{y2, ■),■■■ Yn{yki ■)) 
converges in distribution on the space Dj^k [0, oo) to a vector {Z{yi, •), Z(y2, ■)■,■■■■, Z{yk, •)) 
of independent D-[i[0,oo)- valued components. 

Remark about mean-zero random walks. We have made no assumption 
on the mean b of the random walk. The results are true also for 6 = 0. How- 
ever, in this case the convergence of Yn does not relate to the hydrodynamic 
limit in the same way because Theorem 2.1 is not the correct limit. The 
relevant hydrodynamic limit takes place on the time scale n'^t and the lim- 
iting evolution is governed by the heat equation. For 6 = 0, Theorem 2.1 is 
completely trivial because u{x,t) = uo{x) and y = x. 

Remark about fractional Brownian motion with H = j. There is, of course, 
a result for Brownian motion that corresponds to the random walk re- 
sult of Theorem 2.2. We state here the equilibrium version. Let A > 0. Let 
{Bi{t):i £l} be a countable collection of independent standard Brownian 
motions on R whose initial locations (and, consequently, the locations at 
any fixed time) are those of a homogeneous, rate A Poisson point process on 
R. Fix yen. Let 

Yx{t) = ^ 1{B^{0) < y, B,{t) >y]-Y. l{Bm > y, Bi{t) < y} 

i i 

be the net current of Brownian particles across the point y during time inter- 
val [0,t]. To have sample paths in Dr,[0, oo), we should replace Yx with the 
right-continuous modification 1^ defined by Y^{t) = Yx{t+). This change 
does not affect finite-dimensional distributions. The covariance of Yx is 

EYxis)Yx{t) = ^{^+Vt-Vt^) fors<t. 
V 2vr 
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Theorem 2.4. As X ^ oo, the process X^^^'^Y^ {■) converges in dis- 
tribution on the space Z)r[0,oo) to fractional Brownian motion Z(-) with 
covariance 

EZ{s)Z{t) = ^L(^ + Vt- Vi- s) fors<t. 

The calculations of this paper can be adapted from the random walk 
situation to the Brownian situation and we omit the explicit proof. 

Before turning to the proofs, we want to compare the independent walks 
with Hammersley's process. 

3. A bound on second-order fluctuations in Hammersley's process. In 

this section we look at Hammersley's process from the same perspective 
from which the previous section studied independent random walks. The 
difference is that now there is genuine interaction among the particles. The 
hydrodynamic equation is a nonlinear Hamilton-Jacobi equation. The char- 
acteristics of the equation can meet and form shocks. Currently, we cannot 
prove the equivalent of Theorem 2.2 for Hammersley's process. We can only 
give bounds on the tails of the second-order fluctuations which suggest that 
if there is a limit, it should be on the n^/^ scale. 

The state of Hammersley's process at time t \s zt = {zt{i):i £ 'Zi). De- 
pending on the preferred interpretation, variable zt{i) E R is the location 
of particle labeled i at time t or the height of the interface over site i. The 
dynamics preserves the ordering zt{i — 1) < zt{i). Particles jump to the left, 
according to this rule. If the state at time t is zt = {zt{i) :i (zZ), then par- 
ticle i has an instantaneous rate zt{i) — zt{i — 1) of jumping, independently 
of all other particles. And when particle i jumps, its new location is cho- 
sen uniformly at random from the interval {zt{i — l),zt{i)). This happens 
independently for all particles i. 

This process can be defined in terms of a special graphical construction 
that utilizes the increasing sequences in a space-time Poisson point process, 
see [1, 15] or [18]. 

The process of increment variables rjt = {'r]t{i) :i G Z) is defined as before 

by 

r]t{i) = zt{i) - zt{i - 1), 

and is also known as the "stick process." The dynamics of r]t operates as 
follows. For each i G Z, at rate equal to this stick-breaking event hap- 
pens: pick u uniformly distributed on [0,?7(i)], and replace the state ry with 
the new state 

f Vii) -u, j = i, 
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In other words, represents the stick configuration after a piece of size 

u has been moved from site i to i + This process can be rigorously defined 
on a certain subspace of the fuh product space [0,oo)^, see [15] for details. 

Next we describe one set of hypotheses under which the hydrodynamic 
limit and the trivial fluctuations (1.2) can be proved. Then we state a bound 
on the size of the second-order fluctuations. The setting is again that of a 
sequence of processes , n € N. 

Assumption D. Assume given a nondecreasing Lipschitz function uq 
on R and a bounded, continuous, nonnegative function vq on R. For all 
n, Zq{0) = and the initial increment variables {7/o(*) £ Z} are mutually 
independent with means and variances 

E[r]J^{i)]=nuo(^^^ -nuo(^^—^^ and \aT[r]l^{i)] = vo(^^y 

Furthermore, assume uniformly bounded fourth moments: 

sup E['i]Q{i)^] < oo. 

neN,iGZ 

Let u{x,t) be the unique viscosity solution of the Hamilton-Jacobi equa- 
tion 

(3.1) ut + {uxf = 0, u{x,0) =uo{x). 

Equivalently, u is defined for t > by the Hopf-Lax formula 

{x-yf 



(3.2) uix,t)= inf \uo{y)+ 

y:y<xl 4t 

The hypotheses guarantee that there exists a nonempty compact set I{x, t) C 
(—00,2:] on which the infimum in (3.2) is achieved: 

f (x — y)'^' 
I{x, t) = iy<x: u{x, t) = uo{y) H — — 

A point {x,t) is a shock if I{x,t) is not a singleton. Equivalently, u{x,t) 
is not differentiable in the x variable at (x,t). For a fixed t > 0, there are, 
at most, countably many shocks. Shocks cannot happen for the linear equa- 
tion (2.14) of independent particles, because its characteristics are parallel 
straight lines. 

Here is the starting point: the hydrodynamic limit and the fluctuations 
transported by the characteristics. 

Theorem 3.1. Assume Assumption D. Then, for each {x,t) G R x 
[0,00), 

(3.3) lira z^^{[nx]) = u{x,t) in probability 
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and 

/3 4^ lim f ^nti[nx])-nu{x,t) .^^ zj} {[ny]) - nuojy) \ _^ 



in probability. 

It is possible to construct the processes zf on the same probability space 
with a two-sided standard Brownian motion B(-) such that these limits hold 
in probability, for all {x,t): 

rnrr^ 1- zlit{[nx]) -nu{x,t) fy 

(3.5) hm ^^^^^^ — '-^ ^-^-L= mf B{ vo(s)ds]. 



This theorem is proved in [18]. The hypotheses for (3.5) in [18] are more 
stringent than the ones used above (in [18] the initial increments are assumed 
exponentiahy distributed) and the conclusion is stronger (a.s. convergence). 
The argument in [18] gives convergence in probability in (3.5) under the 
fourth moment bound included in Assumption D. 

The infimum in (3.4) and (3.5) is in some sense the same infimum that 
appears in the Hopf-Lax formula (3.2), which is inherited by a particle-level 
variational formulation (5.2). 

The result of Baik, Deift and Johansson [2] gives the fluctuations for Ham- 
mersley's process from the following particular deterministic initial state: 
Zo(i) = for i < and Zo(i) = oo for i > 0. In this situation the number 
of particles in space interval (0,a;] at time t equals the maximal number 
L(x,t) of space-time Poisson points on an increasing path in the rectangle 
(0, x] X (0, t] . This connection comes from the graphical construction of Ham- 
mersley's process. The distributional limit for n~^/'^{L(nx, nt) — 2n\^xi} in 
[2] can be translated into a limit for n~^^'^ {zntUnx]) — ?i3;^/(4t)} for x,t>0. 

We saw for the independent random walk model that the scale of the fluc- 
tuations from deterministic initial conditions is the scale of the second-order 
fluctuations. So, given the Baik-Deift- Johansson result, we would expect 
the next order of fluctuations for Hammersley's process at scale n^/^. To 
capture these fluctuations, fix {x,t) G R x (0,oo) and define 

yn = {znt{[nx])-nu{x,t)} - inf {z^^ {[ny]) - nuo{y)} . 

yel(x,t) 

We shall prove a bound on the tails of Yn that suggests n^^^ as the correct 
order. We need one more assumption. Given {x,t), let 

Hy) = uo{y)+^-^^ 
be the quantity minimized over y in the Hopf-Lax formula (3.2). 
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Assumption E. Given {x,t), the minimizers in (3.2) are uniformly 
quadratic, in other words, there exist ci,5 > such that 

(3.6) ^y) - <^{y) > ci{y - yf 

for all y G R and y G I{x, t) such that \y — y\ < 6. 

Theorem 3.2. Let {x,t) G R x (0,oo). Under Assumptions D and E, 
the sequence 

Ln^/'^logn J 

is tight. 

We turn to proofs, beginning with the random walk model. 

4. Proofs for the random walk model. The main work is in proving Theo- 
rem 2.2. Along the way we derive an estimate that takes care of Theorem 2.3. 
Last we explain how Theorem 2.1 follows. We start by proving the conver- 
gence of finite-dimensional distributions to the correct limit, and then prove 
tightness at the process level. 

Using (2.20), write Yn{t) = Yn,i{t) - Yn,2it), where 

Yn,i{t) = ^ (0) > [ny] + 1, X^int) < [ny] + [nbt]} and 

(4.1) 

Yn,2it)=J2l{X^m < [ny],Xrint) > [ny] + [nbt]}. 

i 

Yn^i{t) and 1^,2 (i) represent contributions of slow and fast random walks, 
respectively. Next write 

Y„{t) = {EYn,i{t) - EYn,2{t)) 

(4.2) 

+ {Yn,l{t) - EYn,l(t)) - (y„,2(t) - EYn,2{t)). 

We look at the behavior of these three terms on the scale n^/^. Note that 
Yn^i{t) and 1^,2(0 are independent of each other. 

4.1. Convergence of finite- dimensional distributions. 

Proposition 4.1. Fix N time points 

0<ti<t2<-- - <tN- 

As n ^ oo the vector n~^^^{Yn{ti),Yn{t2), ... ,Yn{ti\f)) converges in distri- 
bution to the mean-zero Gaussian random vector {Z{ti), Z{t2), ■ . ■ , Z{tN)) 
with the covariance defined in (2.21). 
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Proof. Since 1^(0) is identically zero, we may as well assume that 
ti > 0. By the Cramer-Wold device, it suffices to show the convergence of 
the linear combination 

N 

(4.3) n-'/^Y.^^^niU) 

1=1 

for an arbitrary vector = {Oi, . . . ,9j\[) ^ and arbitrary G N. 

In Lemma 4.6 below we show that n-'^/^EYn{t) = n-'^/^{EYn^i{t) - EYn,2{t)) 
vanishes as n — > oo. Using the decomposition (4.2) and ignoring the first 
term that vanishes, what we actually prove is the weak convergence of the 
difference 

N N 

(4.4) n~^/^Y.^i{Yn,i{ti) - EYn,i{ti)) - n-^'''Y.e,{Yn,2{U) - EYn,2{ti)). 

1=1 i=l 

Since the two sums above are independent, we can treat them separately. 
We shall show below that they converge to mean-zero normal distributions 
with variances 

N 



E^'V^ Poiy) / PiBu > z)P{Bt, < z)dz 

i=l ^ -^-"^ 

+ vo{y) f P{Bt^<zfdz 



0-1 



(4.5) +2 E 0,ej^2\po{y) f [P{Bu>z)P{Bt^<z) 



-P{Bt^>z>Bt^)]dz 
+ vo{y) f P{Bt^ < z)P{Bt^ < z) dz\ 

J — oo J 



and 



^ f PCX) 

4 = E ^'^^ Poiv) / P{Bu > z)P{Bu < z) dz 
1^1 ^ -^0 

+ vo{y) P{Bt,>zfdz^ 
(4.6) + 2 E ^.%v^|po(y) r[P{Bu > z)P{Bt^ < z) 



l<i<3<N 



-P{Bt^>z>Bt^)]dz 
+ vo{y) P{Bt, > z)P{Bt^ > z) dzj. 
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Bt above stands for standard one-dimensional Brownian motion. The quan- 
tities cjf and (jI are actually equal. Their different formulas connect to the 
calculations by which they arise in the proof. The limit of (4.3) is then a 
centered Gaussian with variance cr^ + al . To evaluate this sum, note that 
the terms with po{y) add up to integrals over (—00,00), while the terms 
with vo{y) are actually equal. We leave the proof of the next lemma to the 
reader (calculus and Fubini's theorem are needed). 

Lemma 4.1. For < s <t, we have these formulas: 



/oo 
P{Bs>z)P{Bt<z)dz 
-00 

(4.8) / P{Bs> z)P{Bt> z)dz 

Jo 



00 V 2'/r 



and 

(4.9) / P(Bs> z>Bt)dz 



2V27r 
^/t^s 



2tx 

From this lemma and the definitions (4.5)-(4.6), we get 
i=i ^ 



vr 



= 0^0,EZ{U)Z{t,), 

where the last equality comes from (2.21). Thus, the linear combination in 
(4.3) converges in distribution to the linear combination X^i^i ^i-^(^j)- Since 
the vector 9 was arbitrary. Proposition 4.1 follows. 

Equations (4.7)-(4.9) can be manipulated to show that the functions 
\/s + t — ^/s V t — s At and ^/s + \/t — y/s + t are positive definite. This en- 
sures that (2.21) is a legitimate covariance of a Gaussian process for all 
nonnegative values po{y) and vo{y). 

It remains now to prove the weak convergence of the sums in (4.4) and 
the vanishing of the mean n~ -^/"^EYnit) in the limit. 

Lemma 4.2. n ^ 00, n~^/'^J2iLi(^i(Xn,i{ti) — EYnA(ti)) converges 
weakly to a mean-zero normal distribution with variance defined by (4.5). 
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Lemma 4.2 will be proved after some preliminary steps. Relabel the ran- 
dom walks so that Xm,j{-) is the jth random walk that starts at site [ny] + m. 
Then 

Yn,iit) = E E HXmjint) - X^j(O) < [nbt] - m}. 

m=l j=l 

Since the random walks are independent of the initial occupation numbers 

oo 

(4.10) EYn,i{t) = E Po (N] + ^) ■ P{X{nt) < [nbt] - m}, 

m=l 

where X{-) represents a random walk with rates p{x) starting at the origin. 
Write 



n - '{Yn,l{t) - EYn,lit)) = E Um{t), 



-1/4 

m=l 

with mean zero summands 



Um{t) = H^mjint) - X„j(0) < [nbt] - m} 

- n-i/^p[}([ny] + m) • P{X{nt) < [nbt] - m}. 
For fixed n and t, the variables {Um{t) :m > 1} are independent. Abbreviate 

N 

Um = GjUmitj)- 
i=l 

Rearrange as follows: 

N oo N oo 

n-^'^Y.^,{Yn,i{ti)-EYn,i{ti))=Y,Y.^,Um{ti)=Y.U^ = Si + S2, 

i=l m=l 1=1 m=l 

where 

[r{n)y/n.] oo 

5*1= E ^™ ^^"^ "^2= E 

"1=1 m=[r(n)v^]+l 

We shall apply the Lindeberg-Feller theorem to Si and show that 5*2 —> 
in L^. To this end, we make r(n) y oo sufficiently slowly. Let 

H{M)= sup E[r,^{xfl{ij^{x)>M}]. 

n>l,x&Z 
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H{M) — > as M — > oo by assumption (2.8). Now choose r(n) = o(\/logn) 
so that r{n) y oo while 

(4.11) r{n)H{n^/^)^0. 

The condition r(n) = o(\/logn ) is imposed so that later we can use assump- 
tion (2.9). _ 

We first show ES2 — > 0. S2 is a sum of independent mean-zero terms Urm 
and so 

00 W 00 

m=[r{n)y^]+l m=[r{n)y/n]+l 

We wrote \\6\\ for the Euclidean norm and used the Schwarz inequality. Since 
N is fixed it suffices to show that for a fixed t, 

00 

(4.12) lim J2 EUl{t) = 0. 

Recall that the variance of a random sum Tk = J2iLi with i.i.d. sum- 
mands Zj independent of the random K is 

(4. 13) Var [Tk] = EK ■ Var Zi + {EZif ■ Var K. 

The variance of the indicator l{X{nt) < [nbt] — m} is P{X{nt) < [nbt] — 
m}P{X{nt) > [nbt] —m}. X{nt) is a sum of a Poisson(nt) distributed num- 
ber of independent jumps, each jump distributed according to {p{x)}. Hence, 
the variance of X{nt) is ntK2. By Donsker's invariance principle, the process 
{{X{nt) — [nbt])/ y^nK2 : t > 0} converges weakly to standard one-dimensional 
Brownian motion. 

By the definition of Um{t) and (4.13), 



E[Uiit)]=Ys.r 



„-i/4 ^ l{X^,j{nt) - XmjiO) < [nbt] - m} 



(4.14) = n-^/'^p^{[ny] + m)P{X{nt) < [nbt] - m}P{X{nt) > [nbt] - m} 
+ n~^^'^Vo{[ny] +m)P{X{nt) < [nbt] -m}^. 
In particular, we get the bound 

(4.15) E[Ul{t)] < Cn~^/'^P{X{nt) < [nbt] - m} 

by the uniform bound (2.8) on the moments. 

By standard large deviation theory and assumption (2.6), as s ^ 00, the 
random walk X(s) has a rate function I which is convex and quadratic 
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around its unique minimum at b. Consequently, for arbitrarily small a > 0, 
there exists a constant < i^' < co such that 

v2 



Kz^, \z\ < a, 
K\z\, \z\ > a. 



(4.16) /(6 + 2)>| 

The rate function gives the bounds 

P{X{s) <sb- su} < exp{-s/(6 - u)} 

and 

P{X{s) >sb + su} < exp{-s/(6 + u)} 

for all u > and s > 0. Even though a large deviations rate function is an 
asymptotic notion, these bounds are valid already for finite s by virtue of 
superadditivity. We use these first in the form 



P{X{nt) < [nbt] -m} < exp|-nt/^6 - 



ni) 

(4.17) 

^ I exp< —K — - >, <m < nta, 
[ exp{— ETm}, m > nta. 

Consequently, applying (4.15) and (4.17), 

oo 

m=[r{n)y^.]+l 

[nta] oo 

<C7n-^/2 exp{-Km^int)-^} + Cn-^^'^ E exp{-Km} 

m=[r{n)^/n]+l m=[nta]+l 
oo , „ C 



< C / e-^^*''--' dx + -r^^ ■ n-V2g-Xant_ 

J r(n.) 



1-e-^ 



This vanishes as n — > oo due to r{n) — > oo. We have proved (4.12). 

Next follows the application of Lindeberg-Feller to S\ . Two things need 
to be checked, namely, that 

(4.18) Ju^ E E^l^{^m\>e]\=^ 

m=l 

for any e > 0, and, second, that 

(4.19) Jim ^Wl\=ol 

m=l 
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(See Theorem (4.5) on page 116 in [5].) Note that 

and so by the uniform bound on moments, 

\Um\<Cn-'^\ri^{[ny]+m) + l). 

(4.18) follows from this and property (4.11) of r(n). In (4.11) we used n^^ 
simply because for any e > 0, en^^'^ > n^^^ for large enough n. 
We turn to verify (4.19). 

TV 

i=l l<i<j<N 

and so the sum in (4.19) can be expressed as 

[r(n)Vn] N 

(4.20) Epl]=Y,efSi,,iU) + 2 ^.^.•5i,2(t.,i,), 

m=l i=l l<i<j<N 

where we abbreviated 

[r(n)yn] 

5i,i(t)= Y EPiit)] 

m=l 



and 



[r(n)Vn] 

Sl,2is,t)= Y E[Um{s)Um{t)]. 
m=l 



Lemma 4.3. Fort>0, 



lim 5i,i(t) = /9o(y) r P{Bt < z)P{Bt > z) dz 
+ voiy)^ P{Bt<zfdz, 



-oo 

where Bt is standard one- dimensional Brownian motion. 
Proof. By (4.14), 

Y E[Ul{t)]=n-^l^ J2 pU[ny]+m)P{Xint)<[nbt]-m} 

m=l m=l 

X P{X{nt) > [nbt] - m] 
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(4.21) 



n 



-1/2 



[r(n)v^] 



J2 ^0 (N] + m)P{X{nt) < [nbt] - m}^ 



m=l 

The last equality above defines the sums Ti^i and Ti^2- We work with Ti^i, 
and leave the analogous arguments for Ti^2 to the reader. 

We bound Ti^i from below. In the calculation that follows, L = L{n) is 
the integer that appeared in assumption (2.9). The o(n~^/^) error term 
below that comes from that assumption is uniform over k because kL = 
0{r{n)^/n) = 0{yjn logn ), which is permitted in assumption (2.9), 

[r(n)Vn] 

Ti,i = p'^iiny] + m)P{X{nt) < [nbt] - m}P{X{nt) > [nbt] - m} 

m=l 

[L-ir(n)v^]-l L 
k=0 j=l 

X P{X{nt) < [nbt] -{k + l)L}P{X{nt) > [nbt] - kL} 

[L-ir(n)v^]-l 

>n-V2 J2 {Lp,{y) + L-o{n-^'^)) 

X P{X{nt) < [nbt] -{k + l)L}P{X{nt) > [nbt] - kL} 

[L-V(n)v^]-l 

> po{y)n^^^'^L J2 P{Xint) < [nbt] - {k + 1)L} 

k=0 

X P{X{nt) > [nbt] - kL} + 0(72"^/^) • 0(r(n)) 



k=0 



+ o{n-^l^)-0{r{n)). 



The last term o(n~^/^) • 0(r(n)) — > as n ^ 00 because r(n) = o{\/\ogn). 
As n — 00, a Riemann sum argument, together with the large deviation 
bounds (4.17), shows that the main part of the lower bound converges to 



Po(y) r P{Bt < z)P{Bt > z) dz. 
J —00 
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Similarly, one derives an upper bound for Ti i with the same limit. This 
proves the convergence of Ti^i. 

We leave the similar argument for Ti_2 to the reader. This completes the 
proof of Lemma 4.3. □ 

Lemma 4.4. For 0< s <t, 
lim 5i,2(s,t) 

n— >co 

(4.22) = po{y)V^ r {P{B, > z)P{Bt < z) - P{B, > z > Bt)} dz 

J —oo 
rO 

+ VQ{y)y/^2 P{Bs<z)P{Bt<z)dz. 

J —oo 

Proof. Formula (4.13) generalizes in the following way. Assume the 
i.i.d. random variables {Zi} are independent of the random nonnegative 
integer and / and g are bounded measurable functions on the state 
space of the {Zi}. Then 



Gov 



Y^f{Z,),j2g{Zi) 

i=l 1=1 



(4.23) = E 



^ f{Zi) - EK ■ Ef{Zi)\lY.g{Zi) - EK ■ Eg{Z^)\ 
.1=1 J u=i J 



= EK ■ Cov[f{Zi),g{Zi)] +YaiK- Ef{Zi) ■ Eg{Zi). 
Applying this gives 

EUn,{s)Um{t) 

= n~^/^p^{[ny] + 'm){P{X{ns) < [nbs] - m, X{nt) < [nbt] - m} 

- P{X{ns) < [nbs] - m} ■ P{X{nt) < [nbt] - m}) 
+ n"^/^Do ([ni/] + m)P{X{ns) < [nbs] - m} ■ P{X{nt) < [nbt] - m}. 
The probabilities in the first term can be rearranged as follows: 
P{X{ns) < [nbs] - m,X{nt) < [nbt] - m} 

- P{X{7is) < [nbs] - m} • P{X{nt) < [nbt] - m} 
= P{X{ns) > [nbs] - m} • P{X{7it) < [nbt] - m} 
- P{X{ns) > [nbs] - m,X{nt) < [nbt] - m}. 
With m = [z^nK2 ] this last expression converges to 

P{Bs > z)P{Bt <z)- P{Bs >z>Bt), 
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which is the integrand of the first integral in (4.22). This points the way, 
and one can fohow the reasoning of the proof of Lemma 4.3. □ 

Together with (4.20), Lemmas 4.3 and 4.4 prove the hmit in (4.19), and, 
thereby, the central limit theorem claimed in Lemma 4.2. Next we need the 
corresponding result for the second sum in the difference (4.4). 

Lemma 4.5. n~^/^YliLi0i{yn,2(ti) — EYn,2{ti)) converges weakly to a 
mean-zero normal distribution with variance o"! defined in (4.6). 

Proof. We have the same argument in principle, 

Yn,2it) =Y1 E ("^) - ^rn,j (0) > [nbt] + m} . 

m=0 j=l 

Then 

oo 

(4.24) EYn,2{t) = ^ p^{[ny] - m) ■ P{X{nt) > [nbt] + m}. 

m=0 

Next write 

oo 

n-^/^(y„,2(t) - EYn,2{t)) = ^rnit), 

m=0 

with independent, mean zero summands 

Vm{t) = HXra,j{nt) - X^j{0) > [nbt] + m} 

i=i 

- n-^^^p'f^{[ny] - m) ■ P{X{nt) > [nbt] + m}. 

With 

N 
i=l 

first separate out the part 

oo 

S2= E ^rn. 

m=[r(n)y/n] 

Use large deviation estimates to show that ^2 ^ in as ?i ^ cx). To the 
remaining part 

[r{n)v'n] — 1 

Sl= Y, ^"i' 
m=0 
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apply Lindeberg-Feller. The details are similar to those in the proof of 
Lemma 4.2. 

□ 

Lemma 4.6. For any <T < oo, 
lim sup n^'^/^\EYn{t)\ = lim sup n~^/'^\EYn,i{t) - EYn,2it)\ = 0. 

Proof. From (4.10) and (4.24), 

EYn,i{t) - EYn,2{t) = S + Ri-R2, 

where 



5 = po(y) < J2 P{^int) < [nbt] - m} - J2 Pi^int) > [nbt] +m}\, 



^ m=l 



m=0 



^1 = E (Poiiny] +m)- po{y))P{X{nt) < [nbt] - m} 



m=l 



and 



^2 = E (Po (N] - "i) - po{y))P{X{nt) > [nbt] + m}. 

m=0 

The part of S in braces equals [nbt] — EX{nt) = [nbt] — nbt. Thus, 15*1 < 
po{y) uniformly over n and t. 

Next we show i?i = o{n}/^) uniformly over t G [0,T]. Ri is a sum of 
bounded terms, and any sum of bounded terms can be rearranged in this 
manner: 



m=l 



L oo 



EE""* 

j=\ m=l 

^ L j-1 ^ L M+j~l ^ L oo 

tEE«- + 7E E 0- + 7E E 



L 



j=l m=l 



L 



j=l m=j 



^ j=lm=M+j 



M 1 L~l 1 ^ 

0W + E1E--+. + 7E E 

m=l j=0 



^ j=l m=M+j 



and so 



00 M ^ L-l 

E "^"^ ~ E ^ E ^rn+j 
m=l m=l j=0 



<0{L)+ E 

m=M+l 
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In our situation L = o(n^/^) from assumption (2.9). We take M = M{n) = 
[cyj n log n ] for a large enough constant c. Then the large deviation esti- 
mates (4.17) show that 

sup 2_\ P{X{nt) < [nbt] — m} — > as n — > oo. 

o<t<T . 

m>cy n \ogn 

It remains to show that the sum 

M , L-1 

^1,1 = E E (/'o (N] +m + j)- pomnXint) < [nbt] - m - j} 

m=l j=0 

is o(n^/^) uniformly over t. Rewrite Ri^i as 

^1.1 =T.lT. (PoiM + m + j)- poiy)) Y.P{Xint) = [nbt] -m-k} 

m=l j=0 k=j 
M oo ^ (L-l)Afc 

= Y.T.P{X(^i) = i^bt]-m-k}- iPoi[ny] + m + j)-pom 

m=l k=0 j=0 
M L-1 ^ k 

= E E^{^(^*) = ["^*]-"^-^}7E(/'o(N]+"i+j)-po(y)) 

m=l fc=0 i=0 
M , L-1 

+ P{X{nt) < [nbt] -m-L}-J2 (Po (Nl + po{y)). 

m=l j=0 

The next to last line above is 0{L) = o(n^/^), as can be seen by replacing 
Po([^y] +"^ + j) ~ Po{y) with a uniform upper bound and then summing the 
probabilities over m. For the last line, use assumption (2.9) to replace each 
^~^J2j'Zo{Poi[''T-y] +'m + j) — po{y)) with o{n~^/'^) uniformly over m. Here 
the assumption m = 0{\/n logn ) is used. Then we have 

AI 

o(""^^^) • E P{Xint) < [nbt] -m-L} = 0(71-^/^) • 0(n^/2) ^ o(^i/4)_ 

m=l 

The last 0(?i^/2) bound comes from 

00 

J2 P{X{nt) < [nbt] -m} = E[{[nbt] - X{nt))^] <Ci + C^n^''^ 

uniformly over t > and n. This is a consequence of the central limit theorem 
and uniform integr ability from assumption (2.6). 

We leave the similar treatment of i?2 to the reader. □ 

Proposition 4.1 is now proved, and we turn to tightness at the process 
level. □ 
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4.2. Moment bound for time increment. Using representation (2.20), the 
difference Yn(t) — Yn{s) for < s < t simplifies to 

Yn{t) - K„(s) = J2{l{X^{ns) > [ny] + [nbs^X^nt) < [ny] + [nbt]} 

i 

- (ns) < [ny] + [nbs],X^{nt) > [ny] + [nbt]}). 

Let 

Yn{t)=Yn{t)-EYn{t). 

By virtue of Lemma 4.6, weak convergence of the processes n~^^'^Yn{-) is 
equivalent to weak convergence of n~^^^Yn{-). We shall, in fact, work with 
the centered processes Yn{-). 

Relabel the random walks so that 1 <j < r]Q{m), are the random 

walks that start at site m in the process indexed by n. Let 

A^j = {X^jins) > [ny] + [nbs],Xli,j{nt) < [ny] + [nbt]}, 

Bm,j = {Xl^jins) < [ny] + [nbs],Xli,j{nt) > [ny] + [nbt]} 

and 

Gm= (1a„„-1b„.,)-P^M(P(A™,i)-P(S^,i)). 

Then 

Yn{t)-Yn{s)= Y Gm 
meZ 

is a sum of independent mean-zero random variables. Recall now assumption 
(2.8), according to which E[j]q{x)^] is uniformly bounded over n and x. 

Proposition 4.2. There exists a constant C such that for all n £N 
and < s <t, 

(4.25) E[{Yn{t) - Yn{s)f] < C(nV2(t _ s)l/2 + ^3/2(^ _ ^^3/2 ^ 

Proof. To prove Proposition 4.2, start with 
E[{Y^{t) -Yr,{s)f] 

(4.26) = E[Gra,G,n,---Grn,] 

mi,...,m6eZ 

= YE[Gl] + (l] Y E[GijE[Gl^] 
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m.i<m2 



6 

2 2 2 



E E[Gl^]E[Gl;,E[Gl^] 

<m2<in3 

<g\y. E[Gl] + ( E EgA ( E EgA 

ImeZ VmeZ / VmeZ / 

+ (t.\egu)\(t.egiX]. 

VmeZ / VmeZ / J 

Above we collected terms E[GmiGm2 " ' " Gme] according to how many times 
distinct sites appear among the indices mi , . . . , mg . Independence and EGm = 

eliminate all terms E[GmiGm2 ■ ■ 'Gmg] where an index appears by itself. 
This point is actually critical for obtaining (4.25). 

Lemma 4.7. There exists a constant G such that for each positive integer 

1 <k <6 and for all m, 

E[\Grn.\''] < G{P{Ara,l) + P{B^,l)). 



Proof. 



E 

1=0 



E 



E iiA 



- 1b 



X p^{7n)''-'\P{Am,l) - P{Bm,lt-\ 

The terms with i < k are bounded by C{P{Am,i) + P(i?m.,i)) as the con- 
clusion demands. For the i = k term we bound as follows: since \P{Amj) — 

PiBm,j)\<l, 



E 



E iiA 

- i=i 



- 1b ■ 



En'?oM = M-i^ 



h=l 

oo 



VIIa -Ib I 

L \j=l 



< E pm^) = h}h^E[\iA^^, - |] 

h=l 
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<E[r]^{m)']{P{A^,,i) + P{Bm,i)). □ 
We can now bound all the sums on line (4.26) by 
(4.27) ^1^"/ ^ ^ E PiAm,i) + C ^ P{Bm,i), 

meZ meZ meZ 

so next we estimate the sums of probabilities on the right. 

Lemma 4.8. There is a finite constant C such that, for all < s <t, 



(4.28) ^\Gm\''<C{^Jn{t-s) + l). 

meZ 

Proof. Write X{-) for a representative random walk that starts at the 
origin. 

Y PiAm,i) = E > N] + [nbs],X!;^^,{nt) < [ny] + [nbt]} 

meZ meZ 

= Y P{X{ns) > [nbs] - m,X{nt) < [nbt] - m}. 

meZ 

Apply the Markov property to turn the sum into 

J2 E Plains) = [nbs] - £}P{X{n{t - s)) < [nbt] - [nbs] +l-m}. 

Rearranging simplifies the sum to 

Y P{X{n{t - s)) < [nbt] - [nbs] + k} 



k<0 



(4.29) 



< Y P{X{n{t - s)) < [nb{t -s)]+k + l} 

k<0 

= E[{X{n{t - s j) - [nb{t - s)] - 1)J 



<CoJn{t-s) + Ci. 



Above we first used [nbt] — [nbs] < [nb{t — s)] + l. Assumption (2.6) gives uni- 
form integrability to all moments of u"^/^(X(ti) — [bu]) as u ^ oo. Due to the 
— 1 inside the expectation, we need the constant Ci in (4.29), CQ^n{t — s) 
alone will not do as t — s — > 0. 

We leave the corresponding calculation for P{Bm,i) to the reader, and 
consider the lemma proved. □ 

We are ready to prove Proposition 4.2. Apply (4.28) to each sum in (4.26), 
remove the squares with 2x^ <x + x^, and let C change its value from line 
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to line: 

E[iYn{t) - Ynis)f] < C{^n{t-s) + 1 + {^n{t-s) + if} 
< C{^n{t-s)+ n^/\t - s)3/2 + 1}. 
Proposition 4.2 is proved. □ 

4.3. Switch to discrete-time process. The processes whose weak conver- 
gence is claimed in Theorem 2.2 are n~^/^l^(t). As observed earlier, it is 
equivalent to prove convergence for the centered processes n~^/^Yn{t). The 
moment estimate in (4.25) is not good enough for tightness, but we can get 
around this by a suitable time discretization. The forthcoming Lemma 4.11 
contains an estimate that gives tightness for the time-discretized process we 
next define. 

Fix two constants a,/3 > such that 

(4.30) f + a</3<|. 

Let 

Wn{t) = n-'/^Yn{n-''[n% 

(4.31) 

= n-i/^(y„(n-^[n'^t]) - EYn{n-^[n'^t])). 

In this section we show that it suffices to prove the weak convergence 
of process Wn by showing that n~^/^Yn and Wn come uniformly close on 
compact time intervals. 



Proposition 4.3. For any < T < oo and e> 0, 

lim pI sup \n^^/^YJt) - WJtM >e\=0. 



Proof. Because the expectations n ^/^EYn{n [n'^f] ) vanish uniformly 
over < t < T by Lemma 4.6, it suffices to prove 



(4.32) 




To prove (4.32), we consider the ways in which Yn{t) can differ from l^(n~^[n^t]) 
some time during [0, T] . First Lemma 4.9 shows that particles that start off 
at least at distance ri^/^+° from [ny] do not contribute to Yn{-) during time 
interval [0,T], in the n — > cxo limit. 
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Lemma 4.9. Let 

Ni{T)= "Y l{X:^j{nt)>[ny] + [nbt] for some 0<t<T} 

(4.33) 

r;(Y{m) 

+ H^mji'^*) ^ ['^y] + [^^*] f(^^ ^^^^ < t < T} 

he the number of particles that start at least at distance n^/^+" from [ny], 
and reach the characteristic some time during [0, nT] . Then for a fixed T, 
ENi{T)^0 as n^QO. 

Proof. We handle the first sum in the definition of Ni{T) and omit the 
similar argument for the other sum. Fix a positive integer M large enough so 
that 1/2 — a(2M — 1) < 0. Again X{-) denotes a random walk starting at the 
origin, and C denotes a constant whose value may change from line to line. 
As Er]Q{m) is uniformly bounded [assumption (2.8)], and by an application 
of Doob's inequality to the martingale X{t) — bt, the expectation of the first 
sum in (4.33) is bounded by 

C V p{ sup {X{nt)-nbt)>e-l\ 

,>„i/2+. l0<t<T J 

<C J2 i'^^' E[{X{nT) - nbT)f^] 

£>„l/2+0! 

£>ni/2+a 

The expectation E[{X{nT) — nbT)"^'^] is 0{n^), as suggested by the central 
limit theorem, due to uniform integrability guaranteed by assumption (2.6). 

□ 

Let 

N2= E 

m: |m-[n3/]|<ni/2+a 

be the number of particles initially within distance n^/'^+°' of [ny]. Fix a 
constant c so that 

lim P{N2 > cni/2+"} =0. 

n— >oo 

The event in (4.32) is contained in the event 

U i sup \Yn{t)-Yn{n~^k)\>ny^e\. 

0<k<[TnP] ^kn-0<t<ik+l)n-0 J 
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For a fixed k, the event in braces implies that at least one of these two 
scenarios takes place: 

(i) At least ^en^^^ particles cross the discretized characteristic s i— > [ny] + 
[bs] during time interval s S [n^~^A;, n^~^(A; + 1)] by jumping. (Note that the 
time interval has been put in the microscopic time scale of the particles.) 
On the event {Ni{T) = 0}, these particles must be among the A''2 particles 
initially within n^'^"*"" of [ny]. Consequently, conditioned on {Ni{T) = 0}, 
the probability of this event is bounded by the probability that A''2 indepen- 
dent rate 1 random walks altogether experience at least ^en^/^ jumps in a 
time interval of length n^~^ . 

(ii) At least ^en^/^ particles cross the discretized characteristic during 
time interval [v}^^k,'n}^^{k + 1)] by staying put while the characteristic 
crosses the location of these particles. For large enough n, the distance be- 
tween the endpoints [ny] + [n^~l^b{k + 1)] and [ny] + [n^~^bk] of the char- 
acteristic is at most 1. Hence, at most 1 site moves from one side of the 
characteristic to the other during this time interval, and so these |en^/^ 
particles must sit on a unique site Xk at time n 

Accounting for all the possibilities gives the bound below. n(cn^/^+"~'^) is 
a mean cj^3/2+a-/3 pgigson random variable and represents the total number 
of jumps among cn^/^"'"" independent particles during a time interval of 
length n^~^, 

p\ sup [Yn(t)-Yn{n-^[n%[>n'/'e\ 
(4.34) < P{Ni {T)>1} + P{N2 > 

+ (P{n(cn3/2+°-/3) > inV^e} + > in^/M)- 

k=0 

The probabilities P{Ni{T) > 1} and P{N2 > cni/2+°} vanish as oo by 
Lemma 4.9 and choice of c. n(cn^/^"*"°~^) is stochastically larger than a sum 

.d. mean 1 Poisson variables, and so a standard large 
deviation estimate gives 

P{n(cn3/2+"-/3) > In^/^e} < exp{-M„/(iM-in^/^e)}, 

where / is the Cramer rate function for the Poisson(l) distribution. By the 
choice of a and /3, M„ > n", while M~^n^^^ oo. Consequently, there are 
constants < Co, Ci < oo, 

[Tn'3] 

^ P{n(cn3/2+"-/3) > in^/^e} < Con^^ exp{-Cin'^) 0. 

To treat the last term in (4.34), we derive a moment estimate for the 
occupation variables uniformly over space and time. 
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Lemma 4.10. Let /c G N and suppose initially sup^gg -^[%("^)'^] < 
Then 

sup E[r]t{m)^] < oo. 

mGZ,t>0 

Proof. Fix x £Z and i > 0, and let 

Cm= H^mjit) =x} 

be the number of particles initially at m who find themselves at x at time 
t. Then 



E[r]t{x)'']=E 



Cm 
^mGZ / 



— -^[CmiCm2 ■ ■ ■ Cmfe] 

r?ii,m2vimfeSZ 

6=1 [mi, 1712,..., mi,) {ki,k2,...,ki,) 

On the last line above we arrange the sum over all /c-tuples (mi, . . . , m^) G Z*^ 
according to the number b of distinct sites among mi, . . . ,mk- The second 
sum on the last line is over 6-tuples (mi,...,m;,) of distinct sites from Z. 
The third sum is over 6-tuples {ki, . . . ,ki,) of positive integers such that 
ki + ■ ■ ■ + kh = k, and (fc^^j^-'-fc;,) counts the number of ways ki mi's, /c2 m2's, 
and so on can be arranged into a /c-tuple. Since mi, . . . ,mb are distinct, 
Cmi , • • • , Cmfc are independent. 

Calculating as in the proof of Lemma 4.7 gives the bound 

E[Ci]<E[7^oini)'']ptim,x), 

where pt{m,x) =pt{0,x — m) is the translation-invariant transition proba- 
bility of the underlying random walk of the particles. Substituting this back 
above gives the upper bound 

EE E ^ ' 



6=1 (mi, 7712, .••,'";,) {ki,k2,...,ki) 

(4.35) 



kik2 ■■■kb 



Y[ptiO,x-mi) \ ■ lY[E[r]oimi)''' 

1=1 J I i=l 
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Holder's inequality, the moment assumption and ki + ■ ■ ■ + kf) = k give 

i=l i=l 

After this, sum the probabilities pt{0,x — rrii) over each index rrii in (4.35). 
This leaves 




which is a constant that depends on A;. □ 
We turn to the last term of (4.34), 

k=0 

Since sup^^^E[r]f{x)^] < oo by the moment hypothesis (2.8) and Lemma 
4.10, and (3 — 3/2 < 0, the right-hand side vanishes as n — > oo. 

We have shown that the right-hand side of the inequality in (4.34) vanishes 
as n — > oo, and, thereby, proved Proposition 4.3. □ 

4.4. Weak convergence. We first verify tightness of the discrete-time pro- 
cesses. Let 

A = (3- 2/3)7(2/3) E(0,i). 

Lemma 4.11. Fix < T < oo. Then there exists a constant C such that 
for all <ti <t < t2 <T and all n, 

(4.36) E[\Wn{t) - Wn{tl)f\Wn{t2) - Wn{t)f] < C{t2 - ti)^+\ 

Proof. If t2 — ti < n~^, then necessarily either [n^t] = [n^ti] or [n^t2] = 
[n^t] . In either case, the left-hand side of (4.36) vanishes and the inequality 
holds trivially. So we may suppose t2 — ti> n~^. 

By the Schwarz inequality and 2xy < + y^, 

E[\Wn{t) - Wn{tl)?\Wn{t2) " Wn{t)?] 

< E[\Wn{t) - Wn{tl)f]'/^E[\Wn{t2) - Wn(,t)f]'/^ 

< E{\Wn{t) - Wn{tl)f] + E[\Wn{t2) - Wn{t)\% 
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Apply (4.25) multiplied by •nr'^l'^ to both terms, ignoring the constant C in 
the front, to get the upper bound 

n\ J n\ nl^ 

(4.37) 



Since t2 — ti>n ^ 



S 75 = t — II H ^ < Z[t2 — tl), 



and also 

l = (n-/5)V/5<(t2-ti)^/^ 
n 

so the first term in (4.37) is bounded by 

1 / [nft]-[n%] \ ^ ^^^^ _ i^)l/2+l/;3 < 2V2T(/^-l)/{2/3)(i2 _ 

n V / 
Apply similar reasoning to the other terms. Note also that 

(t2-tl)=^/'<rl/2-A(i2-ti)l+A 

and 

Collecting terms gives (4.36). □ 

Propositions _4.1 and 4.3 imply that the finite-dimensional distributions 
of the process Wn converge to those of Z defined in Theorem 2.2. This and 
Lemma 4.11 are the hypotheses needed for Theorem 15.6 in [3]. We conclude 
that the processes Wn converge to the process Z on the space Dii[0,T]. 
Proposition 4.3 then implies that the processes n~^^^Yn converge to Z on 
the space -Dr[0, T]. Since T is arbitrary, the convergence holds on -Dr[0, oo). 
This completes the proof of Theorem 2.2. Instead of Theorem 15.6 in [3], 
one can use Theorem 8.8 on page 139 of [6]. 

4.5. Proof of Theorem 2.3. Let Ni{y,T) denote the random variable de- 
fined by (4.33) to display its dependence on y. On the event 

(4.38) {A^i(yi, T) = Ni{y2,T) = ■ ■ ■ = Ni{yk,T) = 0}, 

the processes Yn{yi, ■), Yn{y2, •)>••• y^niVk, Oi restricted to [0,T], depend on 
disjoint collections of independent random walks if n is large enough. The 
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probability of (4.38) converges to 1 as n — > oo by Lemma 4.9. Consequently, 
the restrictions to [0,r] of the processes Yn{yi, •) become independent in the 
limit. To prove the tightness of the joint process n~^/^(l^(yi, ■),Yn{y2, ■),■■■, YniVk, •)) 
on the space L'j^fe[0,oo), apply Theorem 8.8 from page 139 of [6] to the 
discrete-time process {Wn{yi, •),Wn{y2, •)>•••> ^niVk, each component de- 
fined as in (4.31). The proof of Lemma 4.11 can be adapted to the multi- 
variate case. We omit the details, and consider Theorem 2.3 proved. 

4.6. Proof of Theorem 2.1. Limit (2.16) follows from (2.19) and The- 
orem 2.2. Subsequently, the hydro dynamic limit (2.15) follows from limit 
(2.16) and assumption (2.7). 

To prove (2.17), we take a two-sided Brownian motion B{-), and create a 
bi-infinite sequence of random times 

• • • < Tn,-2 < Tn,-1 < = Tn,0 < Tn,l < r„,2 < • • • 

such that 

(4.39) ao"(x) =j2P0 (-) + V^B{n~'Tr,^,) 



m=l 



can be taken as the initial condition. (If x < 0, the sum actually ranges 
over X + 1 < m < 0.) To achieve this, apply the usual Skorokhod embedding 
(see, e.g., Section 7.6 in [5]) to the independent mean-zero random vari- 
ables {7/0(2;) — Poin) - x GZi} and the two-sided Brownian motion Bn{s) = 
v}/'^B{—). Embed {t/q (^) ~ Po(^) > 0} in the positive half of Bn, the re- 
maining random variables in the negative half of Bn. Then the increments 
{Tn,x — Tn,x-i : X G Z} are independent with means 

(4.40) E{Tn,x - Tn,x-i) = VarK(x)] = vo (^^) , 

and we have the equality in distribution of processes 

{Bn{Tn,cc) - Bn{Tn,x-l) : X G Z} = |r?^(x) - (^^^ : X e z|. 

Note that we have been using (2.11) which is assumed for this part of Theo- 
rem 2.1. Now it is clear that do , defined by (4.39), has the right distribution 
to serve as the initial height function. 

Next we observe the central limit theorem for a^. From (4.39), for y G R, 

— — li{n ln,[ny]) 

(4.41) 

+ V^[ (po (^) - Po{s)) ds + 0{n-^'\ 
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One can show that 

hm n Tn,[ny]= / vo{s)ds in probabihty. 

This follows from Chebyshev's inequality and the moment bound 

E[{Tn,m - Tn,m-if] < CE[{r,^{m) - Er,^{m))^] < Ci. 

As n — > OO, the integral term in (4.41) vanishes by the Holder property (2.10) 
of pq. Consequently, 

^o"(M)-nno(,)^^/P^ X 

in probability. Finally, (2.17) follows from this and (2.16). 

5. Proof for Hammersley's process. The proof of Theorem 3.2 is based 
on the approach and estimates derived in [18]. To save space, we refer to 
that paper for all the groundwork. 

We construct the initial configuration by Skorohod's representation, so 
that 

z^{i) =nuo{i/n) +n^/^B{n-^Tn,i), 

where for each fixed n, {r„_j : i G Z} are the hitting times of appropriate ran- 
dom intervals independent of the two-sided Brownian motion B{-). Section 
8 in [18] discusses this construction. For {x,t) G R x [0,cxd), set 

Cr(x)=n-i/2{^-([nx])-nn(x,t)}. 
In particular, in terms of the Brownian motion at t = 0, 

(5.1) Co (y) = B{n-'Tn,[ny]) + n^^\uo{[ny]/n) - no(y)). 

Hammersley's process has a special graphical construction in terms of 
increasing sequences among rate 1 Poisson points on the space-time plane. 
This representation goes back to [9] and [1]. It can be expressed as follows: 

(5.2) ^r(fc) = .inf>o"(^) + rr(^-^)}, 

where r"'*(m) is the minimal positive h such that the space-time rectangle 
(zQ(i),ZQ{i) + h] X (0,t] contains at least m Poisson points on an increasing 
path. This appears as equation (45) in [18]. But observe that compared to 
[18], in the present paper time arguments have become subscripts and space 
indices have become arguments in parentheses. 
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5.1. Upper bound. Following Section 6.1.4 on page 172 in [18] gives < 
Rn,3 + Co where Co is a constant and 

Rn,3= sup irj '\[nx]-[ny\) — y 

yeiix,t) I 4nt J 

By Lemmas 4.4 and 4.5 in [18], P{Rn,3 > Cn^^^logn) ^ if C is fixed large 
enough. 

5.2. Lower bound. Let e > 0. We shall show that for some constant 
< C < oo, Yn > — Cn^/^logn with probability at least 1 — e for large n. 

Let in be the minimal microscopic minimizer for z"j([nx]) defined by 
(59) in [18]. By Lemma 5.3 in [18], dist{n~^in, I{x,t)) ^ in probability. 
Let > be as in Assumption E. Then for large enough n, there exists a 
random y„ G I{x,t) such that 

(5.3) P{\n-hn-yn\<S}>l-e/i. 

Assume now we are on the event in braces in (5.3). Let [01,02] be a 
compact interval that contains the (5-neighborhood of I{x,t). Following the 
calculation in Section 6.1.2 on page 170 in [18] gives 

\2 



Yr,. > r: 



nt 



{[nx] - Z„) - (l^^^nLj + „l/2{^n(^-l _ ^„(^^)} 



Ant 

+ n{^{n-^in) - ^{Vn)} 
> Rn,i + n^^^CSi^'^in) - CSiVn)} + Cin-\in " ny„)2. 
We used the definition 

f ni/r {nx — i)"^ ] 

Rn 1 = min< F 1 {[nx\ — i) : nai <i< na2 }, 

[ Ant J 

and applied assumption (3.6). By Lemma 4.3 in [18], for large enough n, 
(5.4) P{RnA>-Cn^/^logn}>l-e/8 

if C is fixed large enough. Of the lower bound of K„ given above, it remains 
to handle this part: 

n'/^Coin-'in) - Co (yn)} + cin-^in - ny^f 

(5.5) 

>n^/^{S(?i ^Tn,iJ-B{n ^r„j„j^„])} + cin ^{in - nynf - 02- 
The last constant — C2 accounts for the last terms in (5.1), actually of order 

Let [61,62] be an interval such that jQVo{s)ds G [61 + 1,62 — 1] for all 
y G [01,02]- Let 

w{u) = sup{ I i?(s) — B{r) \ : .s,r £ [61, 62], |s — r| <u} 
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be the modulus of continuity of Brownian motion on the interval [61 , 62] of 
length A = b2 — bi. By Levy's theorem, 

limsup(2ulog(j4/n)) ' w{u) = l a.s. 
Choose 5o G (0, 1/2) so that 

(5.6) P{w{u) < 2(nlog(^/n))^/^ for all u<5o}>l- e/16. 

The following law of large numbers holds: for all large enough n, 



(5.7) p| sup n ^r„j„j^] - / vo{s)ds 



<5o\>l- e/32. 



This can be proved as Lemma 8.1 in [18] is proved. The fourth moment 
assumption enters here, because the hitting times satisfy 

E[{Tn,^-T^^^^lf]<CE[r,^{^)% 

On the event (5.3), both ?/„ and n~^in lie in [01,02]. Then on the event 
in (5.7), both n^^Tn^i^ and n~^T„ lie in [61,62] and so fall within the 
range of the event in (5.6). 

Thus, on the intersection of the events in (5.3), (5.6) and (5.7), the right- 
hand side of the inequality in (5.5) is bounded below by 

(5.8) - 2\Tn,i„ - Tn^[ny„]\^^'^ i^og An)^^"^ + cin~'^ {in- nynf -C3. 

The new constant — C3 includes the earlier constant — C2 from (5.5) and the 
case where 

\Tn,in ~ ^n,[nj/n]l < 

as otherwise this factor should be in the denominator inside the logarithm. 
But since ulog{A/u) ^ as m — > 0, the case of small \Tn^i„ — T^jnyJ I can be 
bounded by a single constant. 

The probabilities of the complements of the events in (5.3), (5.4), (5.6) 
and (5.7) add up to less than e/2. On the intersection of these events, Yn 
is bounded below by — Cn^/^logn, plus the expression in (5.8). It remains 
to show that, with probability at least 1 — e/2, the expression in (5.8) is 
bounded below by —Cn^^^logn for yet another constant C. 

We begin with a simple general fact. Suppose {Xi} are nonnegative, in- 
dependent random variables with bounded means and variances. Let < 
a,b < 00, and C >1 + EXi and Ci > Var for all i. Then, 

( e+k-i ^ 

(5.9) P\ sup ^ Xi>3Ck\ <Ci{a + b)n^^^^. 
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To prove this, note first that the event in braces imphes that 

Xi> Cr?!'^ for some 1 < m < (a + h)n^l'^ . 

Then by Chebyshev's inequahty, 

pi sup X,>Cn^^^\ 

[l<m<(a+6)ni/3 ^„2/3<j<(„,+i)„2/3 J 

<P\ sup {Xi-EXi)>T?lA 

Applying (5.9) to |T„,j„ — (these differences are sums of indepen- 

dent, nonnegative increments), we conclude that on the event (5.3), for large 
enough n with probability at least 1 — e/2, 

\Tn,i^ - Tn,{nyA \ < C'd^n - Vn^/'^). 

From this, the expression in (5.8) is bounded below by 

(5.10) - 2C(|i„ - ny„,|i/2 v n^l^){\ogAnfl'' + cxn~\i,, - ny^f - C3. 

This last expression is no less than — Cn^/^(logn)^/^ for a suitable (new) 
constant < C < 00. This a lower bound for Yn with probability at least 
l-e. 
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